High-spin to low-spin and orbital polarization transitions in multiorbital Mott systems 
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We study the interplay of crystal field splitting and Hund coupling in a two-orbital model which 
captures the essential physics of systems with two electrons or holes in the e g shell. We use single site 
dynamical mean field theory with a recently developed impurity solver which is able to access strong 
couplings and low temperatures. The fillings of the orbitals and the location of phase boundaries 
are computed as a function of Coulomb repulsion, exchange coupling and crystal field splitting. We 
find that the Hund coupling can drive the system into a novel Mott insulating phase with vanishing 
orbital susceptibility. Away from half-filling, the crystal field splitting can induce an orbital selective 
Mott state. 
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The Mott metal-insulator transition plays a fundamen- 
tal role in electronic condensed matter physics [l[ ■ Much 
attention has focused on the one-orbital case, in part be- 
cause of its presumed relevance to high temperature su- 
perconductivity and in part because appropriate the- 
oretical tools for the multiorbital case have until recently 
not been available. In most Mott systems, however, more 
than one orbital is relevant Q and the redistribution of 
electrons among different orbitals leads to new phenom- 
ena such as orbital ordering or "orbital selective" Mott 
transitions. Recent studies of nickelates 0] ^titanates ||, 
cobaltates 0, manganates Q, vanadates [1, [ij E3] and 
ruthenates H, 11 , lj, HI have focused interest on the in- 
terplay between the Mott metal-insulator transition and 
orbital degeneracy. A fundamental question in this field, 
relevant in particular to the issue of lattice distortions 
in strongly correlated materials, is the response of multi- 
orbital systems to a perturbation which breaks the orbital 
degeneracy. In this paper, we show that a two orbital 
model with Hund coupling and crystal field splitting ex- 
hibits two fundamentally different Mott phases, one char- 
acterized by a vanishing orbital susceptibility, and one 
adiabatically connected to the band insulating state. We 
characterize these phases in terms of the atomic ground 
states. 

Multiorbital models are more difficult to study both 
because of the larger number of degrees of freedom, 
and because the physically important exchange and 
pair-hopping terms are not easy to treat by standard 
Hubbard- Stratonovich methods [lj|. Weak coupling ap- 
proaches [12[ have been used to show that exchange and 
pair hopping interactions act to suppress the response to 
a crystal field splitting, and some authors have studied 
the model without exchange and pair hopping terms 
but a reliable extension of these results to physically rel- 
evant Slater-Kanamori interactions and the strong cou- 
pling regime has been lacking. 

Dynamical mean field theory (DMFT) provides a non- 
perturbative and computationally tractable framework to 
study correlation effects and has allowed insights into 



the Mott metal-insulator transition [15(. In its sin- 
gle site version, DMFT ignores the momentum depen- 
dence of the self-energy and reduces the original lat- 
tice problem to the self-consistent solution of a quan- 
tum impurity model given by the Hamiltonian Hqi = 
U\oc + #hyb + -ffbath- For multi-orbital models H\ oc — 
Em f mC,U + Ej,fc,J,m U jklm c]c k c\c m , where to = (i, a) 
denotes both orbital and spin indices, and U^ klm some 
general four-fermion interaction, -f/hyb and -ffbath are 
the impurity-bath mixing and bath Hamiltonians, respec- 
tively. While the DMFT approximation simplifies the 
problem enormously (replacing a 3 + 1 dimensional field 
theory by a quantum impurity model plus a self consis- 
tency condition) , the extra complications associated with 
exchange couplings in multiorbital systems have until re- 
cently prohibited extensive numerical work. Interesting 
progress has been made using a finite temperature exact 
diagonalization technique 0, [l3| , but this approach re- 
quires a truncation of Hbath to a small number of levels. 
In Refs. [3, G3] we have introduced a continuous-time 
impurity solver which can handle the general interactions 
in H\ oc . The method, which is free from systematic er- 
rors, is based on a diagrammatic expansion of the parti- 
tion function in the impurity-bath hybridization -ffhyb- 

Here, we employ this solver to study the physically 
relevant case in which the number of electrons matches 
the number of orbitals. The local Hamiltonian is 

H\ oc = - ^2 + y A(m, CT - 712,(7 ) 

a— 1,2 a a 
a=l,2 a 
(J 

with /.t the chemical potential, A the crystal field split- 
ting, U the intra-orbital and U' the inter-orbital Coulomb 
interaction, and J the coefficient of the Hund coupling. 
We adopt the conventional choice of parameters, U' — 
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FIG. 1: Filling of orbital 1 as a function of U for A/i = 
0.2, 0.6, 1 and several values of J/U. The different curves 
for given A correspond (from bottom to top) to J/U = 0, 
(0.01), (0.02), 0.05, 0.1, 0.15, 0.25, respectively. Open (full) 
symbols correspond to metallic (insulating) solutions. The 
metal-insulator transition is characterized by a jump in filling 
and a coexistence region where both insulating and metallic 
solutions exist. Our data show the region of stability of the 
metallic phase. 



U — 2 J, which follows from symmetry considerations for 
d-orbitals in free space and is also assumed to hold in 
solids. With this choice the Hamiltoman (JTJ) is rotation- 
ally invariant in orbital space and the condition for half- 
filling becomes fi = /i 1/2 = |?7 - §J. In the DMFT self- 
consistency loop we use a semi-circular density of states 
of bandwith it (Bethe lattice). The temperature, unless 
otherwise noted, is T/t = 0.02 and we suppress magnetic 
order by averaging over spin up and down in each orbital. 
No sign problem is encountered in the simulations. 

The main result is shown in Fig. [TJ which for several 
values of A and J/U plots the filling per spin, n\, of or- 
bital 1 as a function of interaction strength. The half 
filling, non-magnetic condition implies ri2 = 1 — n\. In 
the T — ► limit, three phases are found: a metallic phase 
(which may have any value of ri\ between and 0.5), an 
orbitally polarized insulator favored by large A and small 
J, and a Mott insulator (with n\ = 0.5 = n-i) favored 
by large U, small A and large J. If U is increased from 
zero to a small value, the orbital splitting either increases 
(small J/U) or decreases (large J/U), consistent with the 
findings of Ref . [lj§] ■ As interaction strength is further in- 
creased, one of several things may happen: at very small 
J/U, n\ continues to decrease, and the system eventually 
undergoes a transition to an orbitally polarized insulator 
(for large A essentially a band insulator) . For somewhat 
larger J/U, the occupancy ri\, after an initial decrease, 
goes through a minimum and begins to increase. At even 
stronger interactions, one then observes a transition ei- 
ther to an orbitally polarized insulator (where n\ may 
take a range of values) or into a special type of insulator 



Mott insulator 
(spin triplet for J/U=0.25) 




FIG. 2: Phase diagram in the plane of crystal field splitting A 
and intraorbital Coulomb repulsion U for indicated values of 
J/U. For J — the phase boundary is a monotonic function 
of A, whereas for J/U > it peaks near A = y/2J (indicated 
by the dotted lines). The insulating state in the region A < 
\/2 J is characterized by a vanishing orbital susceptibility. 



with n\ = 0.5. 

Figure [2] shows the metal-insulator phase diagram in 
the space of crystal field splitting and Coulomb repul- 
sion for several values of J/U. In the absence of a crystal 
field splitting (A = 0), we observe a metal- insulator tran- 
sition at a strongly J-dependent critical U. This finding 
is consistent with data presented in Ref. (l8| . As A is 
increased, the critical U changes. For J = and fixed 
U , ni decreases until the band is emptied and a metal- 
insulator transition occurs. The monotonic decrease of 
the critical U with A at J = is a special case. For 
J > 0, the first effect of a small A is to stabilize the 
metallic phase. Then, at larger A, a reentrant insulating 
phase occurs. We shall show below that this behavior 
arises from the unusual nature of the insulating state at 
J > and small A, which is characterized at T = by a 
strictly vanishing orbital susceptibility. If A is increased 
at large U , this state makes a transition to an orbitally 
polarized insulator at A k \/2J. We therefore plot in 
Fig. [2] the curves A = \/2J as dotted lines, and sug- 
gest that they correspond to the T = phase boundary 
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FIG. 3: Filling of orbital 1 as a function of A for fixed U 
and indicated values of J/U. Top panel: U/t = 6. Open 
(full) symbols correspond to metallic (insulating) solutions. 
Bottom panel: U/t — 9. Here, all solutions are insulating. 
For crystal field splittings smaller than A c = v2 J (indicated 
by a vertical line) the orbital susceptibility in the T — > limit 
is completely suppressed. Solid lines are for (it — 50, dotted 
lines show results for jit — 12.5, 25 and 100, respectively. 



between two distinct insulating states. 

Figure [3] plots the filling of orbital 1 as a function of 
crystal field splitting for fixed U/t and several values of 
J/U. The leftmost curve in the upper panel shows the 
density variation for J/U = 0.02. At A = 0, the model 
is metallic. The rapid variation of n\ with A reflects the 
large, but finite orbital susceptibility of the metal, which 
for this small value of J is strongly enhanced by U. At 
A/t w 0.325 > \/2J, an apparently first order transition 
occurs to the orbitally polarized insulating state, which 
then evolves smoothly (as A is increased) to the band 
insulator (m — > 0). The two larger J values reveal a 
different behavior. For A < A c » \/2J the insulating 
state is characterized by an orbital occupancy which is 
independent of crystal field splitting. Then, an appar- 
ently discontinuous transition occurs to a metallic state 
with a large orbital susceptibility, which at even larger A 
exhibits a first order transition to the orbitally polarized 
insulating state. The lower panel of Fig. [3] shows the be- 
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FIG. 4: Weight of the different eigenstates of H\ oc for U/t = 
6, J/U = 0.05 and A/t = 0.3, 0.5 and 0.7. The smallest 
crystal field splitting corresponds to an insulating state with 
suppressed orbital susceptibility, the intermediate value to a 
metallic state and the largest splitting to a "band insulator" 
(see Fig. . 



havior for larger U, where the model is always insulating. 
Our data for J = exhibit a rapid variation of ri\ with 
A. The slope is set by the inverse of the Kugel-Khomskii 
superexchange ~ t 2 /U; thermal effects are unimportant 
at fit — 50. For J > and small A, the model is insu- 
lating, with a vanishing orbital susceptibility, then (near 
A = v2J) makes a transition to the orbitally polarized 
insulating phase with a differential susceptibility dn\ / dA 
determined by Kugel-Khomskii physics. Note that the 
transition between the two insulators is sharp only at 
T = 0; for T > a rapid (but smooth) crossover occurs. 

To gain insight into these phenomena, we look at the 
contribution to the partition function from the differ- 
ent eigenstates of the local Hamiltonian. H\ oc has 16 
eigenstates, which we number essentially as in Table II 
of Ref. [17]. For the following discussion it is important to 
note that |6), |7) and |8) are the three spin triplet states 
(with energy U — 3 J — 2/x), while 1 10) and |11) are lin- 
ear combinations of the states |t|, 0) and |0,||) with two 
electrons in one orbital and none in the other. The latter 
two states are coupled by the pair hopping and affected 
by the crystal field splitting. Here, we choose them to 
be eigenstates of i?i oc corresponding to the eigenenergies 
U ± V J 2 + 4A 2 - 2f i: (1 + al)- 1 / 2 ^ U, 0) + a±|0,U», 
a± = ±(v J 2 + 4A 2 =p 2A)/J. In particular, we choose 
1 10) to be the eigenstate with lower energy. 

Figure |4] shows the weights of these states for the three 
phases found at U/t = 6, J/U = 0.05 (see Fig. [3J. 
In the small-A phase, the triplet states are occupied, 
with small excursions into states with occupancy 1 or 
3. We therefore call this phase the triplet Mott insu- 
lator. The triplet states of course have one electron in 
each orbital and gain no energy from orbital polarization 
(the remarkable fact is that this feature is preserved af- 
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FIG. 5: Filling m(ji), n 2 {^) for U/t = 4, J/U = 0.25 and 
A/t = 0.4. Full (open) symbols correspond to insulating 
(metallic) solutions. At half-filling (fj./t = 3.5), the system 
is in a triplet Mott insulating state, for 3.9 < (i/t < 4.6 in 
an orbital selective Mott state, and for /J./t > 4.6 metallic in 
both bands. 



ter coupling to the lattice) . In the metallic phase, a large 
number of states is visited, while in the orbitally polar- 
ized insulator, the dominant local state (whose weight 
increases continuously with A) is a singlet (|10)). The 
triplet states are almost completely suppressed in the or- 
bitally polarized phase. The large-?/ insulator-insulator 
transition exhibits the same features, but without the 
intermediate metallic phase, and is therefore also a tran- 
sition between high and low spin states. Comparison of 
the eigenenergies of the spin triplet states and 1 10) show 
that these levels cross at A = \/2J. Thus, the transition 
from triplet Mott insulator to orbitally polarized insula- 
tor occurs at A c = y2J, consistent with our numerical 
data. We also note that the wave- function of state 1 10) 
depends on the ratio J/ A, leading in the large- A limit 
to rii (A) w (J/4A) 2 . In the low spin phase, the or- 
bital susceptibility has therefore two contributions: one 
originating from Kugel-Khomskii physics and one of or- 
der J 2 /A 3 from H\ oc . The latter explains the roundings 
seen in the right most curve of Fig. [31 

We briefly address the issue of the orbital selective 
Mott transition, which provides a mechanism for local 
moment formation in correlated materials, and has been 



the subject of much recent debate [ill ]. Previous studies 
focused on two-orbital models with different band-widths 
and integer number of electrons. We find that in the 
presence of a crystal field splitting, shifting the chemical 
potential can drive the system into an orbital selective 
Mott state, even if the band-widths are equal. Figure [5] 
shows the filling per spin in both orbitals as a function 
of fi, for U/t = 4, J/U = 0.25 and A/t = 0.4. Doping 
occurs first in one of the bands, leaving the other in a 
Mott state with a magnetic moment. Further change of 
the chemical potential drives the second band metallic. 



In conclusion, we have shown that multiorbital impu- 
rity models with realistic couplings can be efficiently sim- 
ulated with the method of Ref. [l7( . We have presented 
numerical evidence, based on single site DMFT calcula- 
tions, for the existence of two distinct Mott insulating 
phases in a half-filled two-orbital model with Hund cou- 
pling and crystal field splitting. At strong interactions 
and J < \/2A, the system, in the T — » limit, is in 
a phase characterized by a vanishing orbital susceptibil- 
ity, and a spin 1 moment on each site. For J > v2A 
an orbitally polarized insulator is found. The exchange 
terms promote insulating behavior at A = but can 
stabilize a metallic phase at values of A for which the 
non-interacting model is a band-insulator. 

It is interesting to compare our results to recent work 
on the bilayer Hubbard model [ijj l20j| . The model 
which these authors study is equivalent to our model with 
U = U' = J , and A replaced by the interlayer hopping. 
In the low energy sector of this model, only four states 
(essentially our three triplets and the pair hopping state 
1 10)) are relevant, and what these authors describe as the 
Mott insulator to band insulator crossover corresponds 
to our transition (apparently sharp at T = 0) between 
triplet Mott insulator and orbitally polarized insulator. 

The existence of two distinct insulating phases raises 
many interesting questions including the theory of an 
insulator with strictly vanishing orbital susceptibility 
(which should exhibit an orbital gauge symmetry) and 
the nature and properties of the different metal-insulator 
transitions. The physics near the "triple point" remains 
to be studied. Our results away from half-filling suggest 
that lightly doped La2Ni04 is in an orbitally selective 
Mott phase. 

The calculations have been performed on the Hreidar 
beowulf cluster at ETH Zurich, using the ALPS-library 
2 if . We thank M. Troyer for the generous allocation of 
computer time, A. Georges and A. Poteryaev for stimu- 
lating discussions and NSF-DMR-040135 for support. 
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